Nanomechanical Spin-Polarizer 
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Torsional oscillations of a free-standing semiconductor beam are shown to cause spin-dependent 
oscillating potentials that spin-polarize an applied charge current in the presence of intentional or 
disorder scattering potentials. We propose several realizations of mechanical spin generators and 
manipulators based on this piezo-spintronic effect. 
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The field of spintronics comprises the search for novel 
logic and sensing devices that employ the electron spin 
degree of freedom by (excess) spin generation and ma- 
nipulation [H]. In the conventional approach spins are 
injected into normal conductors by ferromagnetic metals 
using an applied electrical bias. An alternative method 
is the spin-pumping by a moving magnetization 0, 0]. 
A net-spin generation is possible without involving fer- 
romagnets at all by making use of the spin-orbit inter- 
action, e.g. by a time-dependent gate acting on a two- 
dimensional electron gas [4|, [5J, or the spin Hall effect 
i. 

Spin-transfer by spin-flip scattering in metal struc- 
tures causes mechanical torques 0, H EJ. Mal'shukov 
et. al. [13] predict that a spin-polarized current can 
induce torsional vibrations in a semiconductor beam by 
strain-induced spin-orbit interaction. The same authors 
speculate about a possible reverse effect, viz. that me- 
chanical motion could induce a spin-polarized current. 

In this Letter, we propose a nano-electro-mechanical 
system (NEMS) that generates spins by the coupling 
to torsional oscillations of a free-standing semiconductor 
bridge/beam/rod that is actuated, e.g., by magnetomo- 
tive [3, 111], electrostatic 12| or piezoelectric [l3j forces. 
Subject to an oscillating strain, the spin-orbit interaction 
in the semiconductor generates a spin splitting which, in 
the presence of a bias, leads to a spin current. In analogy 
with piezo-electricity, in which elastic strain induces free 
charges, this can be termed a piezo-spintronic effect. We 
illustrate the physical principle by a conducting wire in 
the electric quantum limit in which only a single quan- 
tized subband is occupied. Subsequently, we generalize 
the results to the multi-channel case. We also demon- 
strate by numerical simulations that the effect survives 
the disorder that can be expected in real systems and 
discuss the conditions under which it can be observed. 

Let us consider a beam with length L and rectangular 
cross section of width d and thickness a (L 3> d{a)) (see 
Fig[TJ that connects two semi-infinite conducting reser- 
voirs. Results can be easily generalized to axially sym- 
metric rods such as catalytically grown nanowires [141 ] - 
The conducting material is a semiconductor that is grown 



Figure 1: A beam consisting of a semiconductor and insulator 
parts (semiconductor layer of thickness a/2 is on top of insu- 
lator layer of the same thickness) is excited by some external 
source into torsional oscillations. A voltage is applied over 
the device in order to detect the mechanically induced spin 
splitting. 



on top of a dielectric. Without loss of generality we as- 
sume here a structure consisting of a conducting medium 
on top of an insulator both being a/2 thick. The lead- 
ing modification of the conduction band Hamiltonian of 
a semiconductor due to lattice strain reads llal: 
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where to* is the effective mass, Uij are elements of the 
strain tensor, and <7j are Pauli matrices. By focusing on 
narrow-gap semiconductors, we may di sreg ard the terms 
proportional to the small parameter (3 |lol |. 

We are interested in the lowest energy vibrations of 
the beam that can be described by an isotropic elastic 
continuum model [13]. Elastic flexural (bending) modes 
cause only diagonal uu strains that couple to the elec- 
trons only via the small (3 term (see Eq. |TJ)) [l(J. The 
strain due to torsional (twisting) vibrations is given by 
El: 



= r(y,t) 



dx 
dx ' 
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where r(y,t) = dip/dy is the derivative of the torsion 
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angle tp with respect to y. The function x characterizes 
the cross-section geometry of the beam and depends here 
only on x and z. It satisfies the equation A\ = — 1 with 
vanishing boundary conditions [17j. We adopt the thin 
plate geometry a <C d which leads to x(x,z) « — (x 2 — 
a 2 /4)/2 and the Hamiltonian: 
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Eq. <(3j) is similar to a Rashba spin-orbit Hamiltonian; 
however, electrons can move here in three dimensions and 
the coupling strength is time and position dependent. 

We now turn to the lowest electronic subband limit, 
disregarding intrinsic spin-orbit interaction, e.g. Rashba 
type, and assuming that the strain induced perturba- 
tion is weak. The free-electron lowest energy states 
read ^(x,y, z) = Ro(x,z)$(y), where Rq(x,z) ~ 
sm(w z/c£) sm(2irx / a) is the lowest subband and $(y) is a 
spinor function. The projected one-dimensional Hamil- 
tonian then reads: 
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where V{y) is the potential due to impurities and we 
disregarded terms ^ (71") 2 . Electrons with up and down 
spins turn out to be uncoupled and subject to effective 

vector potentials of opposite sign, A = 



±h rr{v ' t)a y. 



Since VxA=0, this vector potential does not describe 
an effective magnetic but a spin-dependent electric field: 
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The equation of motion for the torsional angle ip (y, t) of 
the beam reads 



dy 2 9 dt 2 



0, 



(6) 



where I = J(z 2 + x 2 )dzdx ~ ad 3 /12 is the moment 
of inertia of the cross-section about its center of mass, 
p the mass density and C is an elastic constant de- 
fined by the shape and material of the cantilever. C — 
^pda 3 for a plate with a <§; d, and p, is the Lame 
constant. The general solution of Eq. iJBJ is a plane 
wave tp = ipoe lu>t±lky , where k = lu/c is the wave num- 
ber, c = 2cta/d = y/C/(pI), and ct = yjp.jp is the 
sound velocity. Throughout this paper, we consider 
a doubly-clamped beam in which the lowest harmonic 
tp = ipa sin(fey) sin(wt) is excited, where u> — ck and 
k = tt/L is the wave number (see Fig. Q]). The stand- 
ing mechanical wave creates an oscillating electric field 
E = dA/dt that is exactly out of phase for spin-up 
and spin-down electrons. In the Born-Oppenheimer ap- 
proximation, the strain induces a parametric potential 
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Figure 2: Maximum spin-polarization in units of the effective 
splitting Uo/Ef as a function of the width of the electron 
waveguide in Fig. [T] The frequency of the mechanical oscil- 
lations uj — 10GHz which corresponds to Uo/Ef = 2 x 10~ 4 
(length of the rod L = 1/im, Fermi length = 30nm, the 
delta function strength v 2 m* /ft 2 = OAEf and m* = 0.06m). 

U(y) = JZf^g9 sin(yy) cos(ujt) (cf. Eq. ©) that is 

adiabatically followed by the electrons. 

The strain-induced potentials ±U(y) vary only slowly 
and do not yet spin-polarize a charge current signif- 
icantly. However, defect scattering strongly amplifies 
the piezo-spintronic effect, as illustrated now by a sin- 
gle short-range potential scatterer V{y) = v8(y — L/2) 
located in the middle of the beam. Disregarding the small 
intrinsic effect caused by ±U(y), the probability that an 
electron with Fermi wave number hp is transmitted by 
the scatterer reads: 
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where Uq = . According to the Landauer con- 
ductance formula, the spin polarization of a charge cur- 
rent becomes: 



P = 



r t -r 
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(8) 

where the Fermi energy Ep — h 2 k 2 F /2m* . The spin- 
polarization oscillates in time with the beam vibration 
frequency. 

Generalization to a multichannel wire is facilitated by 
the following gauge transformation 
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with e i/ ^' y,t ) cr - s = tcos(/) + id z sin(f), which leads to 
the transformed Hamiltonian: 
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It is convenient to choose 
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allow many occupied subbands along the z axis but re- 
strict considerations to the lowest subband along the x 
axis, which is the case for a laterally weakly confined two- 
dimensional electron gas. After projecting Eq. (flO|) to 
the lowest mode in the x direction, we obtain the follow- 
ing two-dimensional Hamiltonian: 
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can be 
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interpreted as a spin precession length and V(y, z) de- 
scribes two-dimensional disorder scattering. r(y, t) is 
here still arbitrary, but we limit our attention to the low- 
est vibrational mode as before. 

The terms proportional to cr x ( y )k z in Eq. lfl~2]l induce 
subband transitions; however, these do not affect trans- 
port when the precession length l so is larger than the 
width of the channel. In the limit of a long and nar- 
row beam, we may again treat the time dependence of 
the Hamiltonian Eq. (fl2|) parametrically in terms of 

d hwL 
the frequency uj. In the limit A so - — <C Uq = 



we can further simplify Eq. (f!2| by disregarding sub- 
band transitions. A simple short-range wall potential 
V(y,z) — v5(y — L/2) does not lead to subband transi- 
tions either. Our system then reduces to a collection of 
independent channels, which lead to a total spin-current 
polarization 



M 



M 



P = ^(T mT -T m; )/^(T mT +T mi ), (13) 

m m 

where m is the index and M the total number of trans- 
port channels. Here, 



2(kl/2m*±U co S (mt)) 
mT(i) 2{h 2 k 2 m /2m* ± U Q cob(w*)) + v 2 m*/h 2 



, (14) 



where k m is the wave number of an electron in the channel 
m at the Fermi energy. In Fig. |(2|), we present results 
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Time (in units 

Figure 3: Spin-current polarization in units of the splitting 
parameter Uo /Ef as a function of time for the torsional oscil- 
lations of the rod in Fig. [TJ The dashed line gives the results 
according to Eq. (1131) . whereas the bold line represent the 
results of numerical simulations based on the tight binding 
model. Parameters are the same as in Fig. [2] 



of Eq. (|13p for the maximum mechanically induced spin- 
current polarization as a function of the beam width. The 
dashed line in Fig. ((3j) are the results of Eq. (fl~3|) for the 
spin-polarization as a function of time. 

The idealized model above allowed us to illustrate the 
physics of piezo-spintronics. We now consider a more re- 
alistic model, including many electron modes, subband 
mixing and arbitrary forms of the potential V(y, z). We 
numerically calculate the scattering matrix using the re- 
cursive Green's function technique and the tight-binding 
representation of the Hamiltonian lfT2]) : 



ij<? C ija C ije 
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where €ij a is the on-site energy that includes V and 
U, t = h 2 /(2m*b 2 ) is the hopping energy and t so — 
h 2 /(2l so m*b) is the hopping energy due to the spin-orbit 
interaction, in terms of the tight-binding lattice spacing 
b. The bold lines in Fig. [3] display our numerical results 
for the polarization as a function of time for the short- 
range wall potential used above. We find good agreement 
with the analytical results for large aspect ratios of the 
beam, as expected. Deviation from the analytical results 

become noticeable when A so ^ — — , 

a 

Finally, we model the potential V(y, z) by on-site An- 
derson disorder with energies distributed over a band 
width W. The beam is represented as a 150 x 9 dis- 
crete lattice and an ensemble averaging over 20, 000 im- 
purity configurations is carried out. A single realiza- 
tion (without averaging) behaves similar to our single 
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Figure 4: Spin-current polarization in units of the effective 
splitting Uq/Ef as a function of time for the torsional oscil- 
lations of the beam in Fig. JT]). Parameters are the same as 
in the figures above. 

defect result in Fig. [3l Averaged results are presented 
in Fig. [H The Anderson disorder strength can be 
measured by an effective 2D mean free path [l8| as 
1-2D = (6\ F E 2 F )/(ir 3 a 2 W 2 ). Here we consider l 2D = 8L; 
3.6L and 2L. When the Anderson disorder is weak, the 
spin polarization is almost a harmonic function of time. 
We may conclude that the piezo-spintronic effect is very 
robust, persisting in a disordered system, confirming the 
qualitative behavior of the analytical model with one 
dominating defect scatterer. 

The choice of parameters above is motivated by the 
following estimates made for a silicon cantilever of size 
a x d x L = (0.05 x 0.15 x l)/im. The resonant fre- 
quency is lu = c t n/L = 10 GHz, using the density 
of silicon p — 2 x 10 3 kgm~ 3 and the Lame constant 
fi = 100 GPa. The maximum angle of torsion (po can 
be estimated by equating the energy dissipation during 
a cycle 2ttlpqC/(QL), where Q is the mechanical qual- 
ity factor, with the energy input rate <PqT, where T is 
the actuating torque applied at the center of the beam. 
Electrostatic torques of T ~ 10~ 12 10~ 15 Nm have been 
already realized [H]. Taking Q = 500, T = 10- 15 Nm 
and C = 10~ 18 Nm 2 , we find <p = 0.2 rad. The corre- 
sponding spin-orbit precession length is l so — lfxm and 

HujL „ 
the spin splitting is Uo = —, — ~ 6 x 10 D eV, using 

the bulk strain-spin-orbit coupling parameter for GaAs 
7 = 2 x 10 8 m _1 [15j which is not expected to be affected 
strongly by the finite structure confinement [16]. This 
leads to polarizations of the order of P ~ 10 -4 . In the 
presence of an applied DC voltage, the mechanical mo- 
tion generates AC spin current. Alternatively, one can 
apply AC voltage synchronized with the mechanical mo- 
tion to obtain a DC spin current. The thus created spin 
accumulation can be detected by e.g. a ferromagnetic 
side contact or by the optical Kerr rotation (2TJ |. 

In order to increase the polarization, one can use semi- 
conductors with lower doping (smaller Fermi energy). 



The width of the rod can be tuned to capture the reso- 
nant features in Fig. [2] Rods with higher quality factors 
can have larger amplitude of oscillations leading to higher 
polarizations. 

Summarizing, we propose a piezo-spintronic effect that 
is based on strain-induced coupling of the electron spin 
degree of freedom and mechanical vibrations in free 
standing semiconductor nanobeams. We show that time- 
dependent strain due to torsional mechanical oscillations 
can lead to a measurable spin polarization of an applied 
charge current. Mechanically generated spin-dependent 
potentials (mechanically-induced Zeeman splittings) can 
be also used for the manipulation of an applied spin cur- 
rents. We propose ways to measure and increase such 
mechanically generated polarization that can be used for 
effective spin injection in spintronic based devices. 
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